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Abstract 

We evaluate the coefficients of the leadmg poles of the complete 
two-loop quark self-energy in the Coulomb gauge. Working in the 
framework of split dimensional regularization, with complex regulat- 
ing parameters a and n/2 — a for the energy and space components of 
the loop momentum, respectively, we find that split dimensional regu- 
larization leads to well-defined two-loop integrals, and that the overall 
coefficient of the leading pole term for is strictly local. Extensive 
tables showing the pole parts of one- and two-loop Coulomb integrals 
are given. We also comment on some general implications of split di- 
mensional regularization, discussing in particular the limit o" — > 1/2 
and the subleading terms in the e-expansion of noncovariant integrals. 



1 Introduction 



Despite serious efforts during the past twenty years to place the Coulomb 
gauge on the same rigorous footing as covariant gauges, we still have no con- 
sistent rules for renormalizing non-Abelian theories in this gauge. There can 
be little doubt, however, that the Coulomb gauge is superior to other gauges 
in at least two respects, namely in the treatment of bound states and in the 
study of confinement in QCD. Since these topics play a crucial role in our 
understanding of the strong interactions, there is clearly a need to put this 
gauge on a sound theoretical basis. To gain a better understanding of the 
various advantages, as well as disadvantages, of the Coulomb gauge, we refer 
the reader to the vast literature on the subject. (For some recent references 
see, for example, refs. 

Unfortunately, progress in the Coulomb gauge continues to be hampered 
by the operator ordering problem in the quantum Hamiltonian, as noted 
by Schwinger in 1962 |T0|. The ordering problem was later re-examined by 



Christ and Lee who demonstrated that the quantum Hamiltonian differed 
from the classical Hamiltonian by special Coulomb interaction terms, la- 
belled (Vi + V2) Inl. A few years later, Cheng and Tsai |jl2| pointed out 
that the (Vi + V2) -terms are equivalent to a distinct class of integrals, called 



energy integrals, which lead to two different types of divergences |T^, |T3[ : 

(a) ordinary UV divergences associated with the structure of space-time, and 

(b) divergences characteristic of the Coulomb gauge, arising from the inte- 
gration over the energy variable go in integrals of the form 

00 r, 

dqi 



It is the divergences from such energy integrals that give rise to ambiguities 
in the Coulomb gauge. These ambiguities come as no surprise, since the 
Coulomb-gauge condition 

V ■ A{x) = (1) 

does not fix the gauge completely, but leaves a residual gauge freedom for 
gauge transformations g{t) that do not depend on the space coordinates x. 
In 1987, Doust and Taylor [ll4| , pTSf came to the conclusion that standard di- 
mensional regularization is incapable of regulating both types of divergences 
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simult aneouslyQ. 

Renormalization in the Coulomb gauge has recently been examined by Baulieu 
and Zwanziger 0] , who treated the Coulomb gauge as the singular limit of the 
Landau- Coulomb interpolating gauge. Zwanziger also exploited the Coulomb 
gauge in the study of confinement [^. Employing a diagrammatic represen- 
tation, he showed that the problematic energy integrals cancel, at least to 
one-loop order. Based on this observation, he assigned to these integrals the 
value zero. However, it is not obvious that this cancellation procedure can 
also be carried out explicitly beyond one-loop order |TH . 



Accordingly, it seems desirable to regulate the Coulomb-gauge integrals on 
an individual basis. To this effect, a novel technique called split dimensional 
regularization was introduced by one of the authors []TB|. The idea is to 
replace the measure d"'q by 

where a and u are understood as parameters in the complex plane satisfying 
a + u = n/2, and the limits a — > 1/2 and u — >■ 3/2 are taken only at the 
end of integration. One may think of split dimensional regularization as a 
special form of dimensional regularization, the special feature being that the 
dimension of the energy component is explicitly specified to be non-integer. 
Whereas in [^] the ambiguous energy integrals were defined to be zero, in the 
context of split dimensional regularization these integrals turn out to be zero 
in a natural way. 

To date, split dimensional regularization has been tested at one loop for the 



gluon self-energy [|16] , and the quark self-energy and quark-gluon three-point 
function [|17|. All Coulomb-gauge integrals appearing in these calculations 
are free of ambiguities and respect the appropriate Ward identities. 
We note in passing that split dimensional regularization has also been applied 
in the context of non-relativistic QCD [Q. 

However, as already alluded to in ref. [IT^, the real challenge comes from 



energy integrals at two loops and beyond [T^. In order to investigate the 



associated ambiguities, we have evaluated the contributions of the leading 
1/e^ poles of the complete two-loop quark self-energy. Working in the frame- 
work of split dimensional regularization, we find that all integrals can be 

^These authors assume that the energy component qo is one-dimensional, while q is 
(n — l)-dimensional. 
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calculated consistently. 

But, whereas the coefficients of the leading poles can be evaluated explicitly, 
the coefficients of the subleading poles are generally more difficult to com- 
pute, since many of the parameter integrals can no longer be expressed in 
closed form. For this reason, the subleading poles have not been evaluated 
in their entirety in this paper. 

It is also worth noting at this stage that the subleading poles exhibit the 
following interesting feature: the regulators a and uo, characterizing the time 
and space components, respectively, appear independently in the various F- 
functions, a clear indication of the special role played by the time component 
in the Coulomb gauge. (See Section |2.3| for more details.) By comparison, in 
the F-functions for the leading poles, the regulators a and uj always appear 
in the combination a + uj = n/2. 

The paper is organized as follows. In Section 2, we first discuss a specific 
integral for which standard dimensional regularization fails to regulate the 
divergence from the energy component, if the latter component is assumed 
to be one-dimensional, and then show how split dimensional regularization 
cures the ambiguity. In the second part of Section 2, we concentrate on two- 
loop integrals, analyzing in particular the properties of the subleading poles. 
Then we comment on some general properties of split dimensional regulariza- 
tion. We demonstrate that the limit cr — > 1/2 a/ter integration is always well 
defined and discuss the implications of split dimensional regularization for 
the subleading terms in the e-expansion of noncovariant integrals. In Section 
3, we outline the calculation of the leading divergence of the two-loop quark 
self-energy and present our results. The highlights of this paper are summa- 
rized in Section 4. Appendix A contains results for the pole parts of one-loop 
integrals (for both integer and non- integer powers of propagators), while Ap- 
pendix B shows the results for the leading poles of several two-loop integrals. 
In Appendix C and D we give general formulas for Coulomb integrals at one 
and two loops in Feynman parameter space. 



2 Split dimensional regularization 
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2.1 Problems with standard dimensional regulariza- 
tion 

The gluon propagator in the Coulomb gauge, given by [IT 



9 

d'i^uil) = 9f,u + -^QtJiQv - ^ {qt^ny + n^q^) , (3) 
= (1,0,0,0), 

is seen to contain the factor While most Coulomb-gauge integrals can 

be computed consistently with standard dimensional regularization, it is no 
secret that the appearance of the noncovariant factor l/g^ in the integrand 
may, under certain circumstances, lead to spurious singularities that require 
special attention. The presence of two or more of these factors^, in combina- 
tion with powers of go in the numerator, is particularly troublesome and can 
lead to the ambiguous energy integrals mentioned in Section 1. To illustrate 
this point, consider the integral 

We shall first show that standard dimensional regularization fails to regulate 
the singularity in the energy component of I^^, (i.e. in /qo), if the energy 
component go of a vector g^ is assumed to be one-dimensional. After Feynman 
parametrization, we obtain 

where 



^Note that in Abelian gauge theories two noncovariant factors containing the same 
loop momentum cannot occur if the number of external legs is < 3, i.e. in self-energy and 
vertex diagrams. 
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Integration over n = (4 — 2e)-dimensional momentum space leads to 

V = tTi^ f' dx f' dyy\DetA\UT{3-'^){VA),{VA)4M^ -VAV]"^-^ 
Jo Jo ^ z 

+h{2-^)A,AM'-VAVr^-'}. (7) 

The important term in Eq. (|^ is DetA. If we assume that the energy compo- 
nent is one-dimensional, while the identity matrix contained in A is defined 
in {n — 1) dimensions, we obtain iDetAja = (x?/)~2, so that 

= irc^ip^)'^ / dxx~2 / dyy2 
Jo Jo 

{r(l + e)(p^ - pon^)iPu - Poriu) - yY~' 

y-'{i-yr}. (8) 

Thus, we see that the "energy component" /qo is ill defined, since it leads to 
the parameter integral 

Jo 

The example above demonstrates that the application of standard dimen- 
sional regularization in the noncovariant Coulomb gauge leads to unavoid- 
able difficulties, as long as the energy component go is assumed to be one- 
dimensional. This conclusion raises another question, namely, how is the 
dimensionality of go defined in standard dimensional regularization? The 
answer to this question is that in general, the dimension of go need not to be 
specified. In covariant integrals, we could distribute arbitrary fractions of e 
to the individual components of a Lorentz vector, as, for example, in 

y"rf"g ^ J c/"«gorf"^girf"^g2rf"~"°""^~"'g3 . 

Due to covariance, the result will invariably depend only on n. However, 
in noncovariant gauges, the situation is quite different. In our example, 
we could as well have assumed that the identity matrix contained in A is 
three-dimensional, thereby leading to an energy component which is (n — 3)- 
dimensional, so that 

|DetA|^ = = (9) 



+5r(.) 



xy 
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In that case, the x-integral in Jqo is given by Jq dx x~2^''^ which is well de- 
fined by invoking analytic continuation in the context of dimensional reg- 
ularization. This is exactly the idea of split dimensional regularization. It 
asserts that in the noncovariant Coulomb gauge, the energy component go of a 
Lorentz vector (as a loop momentum) cannot be treated as one-dimensional. 
Instead, a nonzero "fraction" of the regulator e in n = 4 — 2e has to be 
assigned to the energy component, for instance by defining 

o- = , e^ = c„-e, (10) 

such that a + uj = n/2. The parameter Co- is arbitrary^, as long as we do 
not set it to zero before all parameter integrals have been carried out. After 
integration, we can define the integral at a = 1/2 by analytic continuation, 
invoking the same arguments as in standard dimensional regularization. 
Thus, the technique of spht dimensional regularization imposes the condition 
that the energy component go must be of non-integer dimension. 
Using split dimensional regularization, we find for the integral 1^,^: 

V = iT^'^{f)-'n^n,V{e)-a l^xx'^-^ l\yy-^-'{l-yr (H) 
+ ViVj ~" ^-iid 5ij — terms ; ^, j G {1, 2, 3} , 

where (a + uj) has been replaced byr?,/2 = 2 — e whenever cr, uj occur in this 
particular combination. Contraction of the result (|TlD with n^n'^ yields for 
the pole part of /qq = / d"- q / {q^ q'^ [p — g)^} , 

divJoo = div{-i7r§(p2)-T(e)Beta(l-(T-e,l-e)} (12) 
= i-K - (—2) . 

As a crosscheck, we may write 

div/oo = div / rf^g^l-t^ = div / d^q—^—^ = zvr^- (-2) . 
J q^^ q'^ip — q)"^ J q'^ [p — q)^ t 



■^Note that in [Q, a special version of split dimensional regularization with Ca = 1/2 
has been used. 
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Note that / rf^'^go / -^j^ — Cr^r = 
J J q'^'yp — q)^ 

in the context of spht dimensional regularization, by the same arguments as 
/ dP'q = in standard dimensional regularization. The breaking of covariance 
is manifest in Eq. ([1^) from the function Beta(l — a — e, 1 — e), which contains 
0", rather than a + uj, in the energy component of I^v It is evident from this 
example that split dimensional regularization leads to well-defined parameter 
integrals. 



2.2 Subleading poles 

In order to extract interesting quantities such as the two-loop renormaliza- 
tion constant for the quark self-energy in the Coulomb gauge, it is of course 
necessary to determine not only the leading (1/e^) poles, but also the sub- 
leading (1/e) poles of the two-loop integrals. Since the expressions for some 
of the Coulomb integrals occurring in the two-loop quark self-energy are quite 
involved (as already noted above), the subleading poles of these integrals can 
no longer be extracted in the form of analytic functions. Instead, they have to 
be expressed as combinations of infinite series and hyper geometric functions. 
We shall illustrate the extraction of the subleading poles by means of two 
examples. The first example consists of an integral Ii where all integrations 
can be done in closed form: 

= K I "^^^ ^ Tie) Betafc^ - 1, 1 - e) 

m2 J ki{p-kiY{kiy 

= (p2)-2^r(e)r(2e)-G(e,a,u;), (13) 



r(l - a)T\uj - l)r(l - e)r(l - 2e) 
T{1 - a + e)T{u - e)T{uj - 2e) 



G{e,a,uj) 



where a + u = 2 — e has been used in Eq. (p!3|) . But note that a and u also 
occur separately in G{e, a, uj), not only in the combination a + u. In order to 
extract the leading and subleading poles in e, we use Eq. (p]0| ) and expand^ 

''A factor r^(l + e) has been extracted in order to avoid the Euler je in the expanded 
result. 
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in e. This yields 

h = {f)-^' T\l + e) II + ^ [5 - - 2 log (2)] + finite} . 



(14) 



We will comment on the appearance of the parameter c^, in the sub lead- 
ing pole term, in Section |2]^. Note that the breaking of covariance is also 
manifest in the factor (p^)"^*^. 

As a second example, consider the integral 



(-l)"T(e)Beta(l -e, 1 -e^ ^ 



r-l 



r(e) Beta(l - e, 1 - e)r(-2e) 



^TT 2 J (A;2)i+^(p - k2fki 
r(l + 2e) 



rfi 







duu'^~'2Fi{l + 2t,~2e]l-e-B{u)) ; B{u) 



1 -*2 

Pq — up 

p2 



-p 
1 

e 



p" 



^1 ^ \P 
Using Eq. (0), we obtain 

h = (V)"^-^T2(l + 6){-^ 

1 r °° 1 

e L ^1 \P 



Po 



+ finite} . (15) 



-1 2^ll-m, -, -, 



2'2 



finite} . (16) 



Again, we see that the breaking of Lorentz covariance manifests itself in 
terms such as p"^ /pi and in the fact that Ca appears in the subleading pole 
term. 

It should be clear from the result for the subleading poles of I2 in Eq. (|T6D 
that a complete determination of the subleading poles in the two-loop quark 
self-energy is beyond the scope of this paper. Afterall, I2 is still a relatively 
simple Coulomb integral since it has only one noncovariant denominator. 
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2.3 Comment on noncovariance in the context of di- 
mensional regular izat ion 

In the context of split dimensional regularization, as already explained in 
Section 2A, the dimension d{qo) of the energy component of an ra-dimensional 
vector is given by d{qo) = 2cr, where (cf. Eq. (p!0|)) 



(17) 



and ea serves to regulate the spurious divergences in the energy component, 
inherent in the Coulomb gauge. The important point is that the limit 
in the energy integrals after integration always exists. This means that 
the "singularities" in the energy integrals, which would be unregulated for 
e^r = 0, are really spurious. They do not show up as poles in e^. Only 
the "usual" UV and IR poles of Coulomb-gauge integrals appear as poles in 
e = 2 — {a + uj), whereas the potential singularities in the energy component 
never show up as r(±/c (1/2 — a)), but rather as r(l — k — a), where k is a 
positive integer. Hence, the limit — *^ with e fixed, i.e. the limit — »■ 0, 
indeed exists. This observation follows from the expressions in Feynman 
parameter space for general Coulomb-gauge integrals given in Appendix 
and 0. First, consider a scalar one-loop integral with r (r > 0) covariant and 
m noncovariant denominators in n = 2{(t + u) dimensions (cf. Eq. (^iD): 



lr+m{0lj, Pi) 



jn r m 

'nn 



2 i\ f^^ (q - Pj)'^''^ (q- Si) 



^ ^ uu r(«.) nz. m) Jo M ' ' [h 

m r m 

1=1 j=i 1=1 

POJ + + ^jPj 

j=l ^j=i -^o j=l 



3=1 1=1 



1=1 
r 771 



;i8) 



A(a,/?) 



j=i 1=1 
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Concerning the appearance of the parameters a and cj, we see that the overall 
F-function (which indicates the UV behaviour) , as well as the term in square 
brackets containing the momentum dependence, always depend on the sum 
(T + a; = n/2. Only the integration over the Feynman parameters Xj (which 
are associated with the covariant denominators) leads to the isolated param- 
eter 0", the reason being that the Feynman parameters yi never multiply a go 
component. 

Referring to Eqs. (^) - (|39D in Appendix we observe that the most dan- 
gerous case for a potential singularity at a = 1/2 occurs when r = 1, / > 2, 
and with factors of Qq in the numerator, that is, in the following type of 
integral (for simplicity, we assume (5i = 1 for all /, since more general powers 
(3i do not spoil the argument): 

/Jra„ m b 

= (-1)"—^/ dxx--'-'^ l[dm6il-x-J2yi) 
i [a) JO J i^-^ i^-^ 

/ dz 2;"+™- 1-1 exp{-2; [xp^ + Y.yi ^1' - (3^^+ E yi ^if]} 

•^0 1=1 1=1 

Carrying out the derivative with respect to 04, we see that the resulting 
x-integral will be of the form 

'rfa;a;—V(x, 2/0, where 7 = ll^^ wo , (20) 







(6 - l)/2 for h odd 



such that 7 is always an integer. The function f{x,yi) corresponds to the 
term in square brackets in Eq. (^) and will, after integration over z, always 
contain the sum a + a; in the exponent. Since f{x,yi) in general does not 
contribute to the leading behaviour for x — * 0, it is irrelevant for our consid- 
erations]^. 

Note that values of a > 1 actually improve the behaviour for x — > in 



^In the special cases where f{x,yi) does contribute to the leading behaviour for a: ^ 
(for example, if / = l,r = l,p = 0), it also provides a regulator e for the exponent of x 
and thus is harmless. 
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Eq. (ppp. By contrast, it is the more severe ultraviolet hehavionr , induced by 
powers of Qq {h > 0) in the numerator, that can render the integral over x 
singular unless the parameter a is complex^. 

Generalization of the above arguments to r > 1 is straightforward. If r is 
the number of covariant propagators, b the power of go in the numerator 
and I > 2, then integration over one of the Feynman parameters Xi requires 
a complex value of a, provided r — 6/2 < is fulfilled (or, respectively, 
I]j=i dj — b/2 < for general propagator powers). 

Summarizing, we conclude from the analysis above that all parameter inte- 
grals are well defined if a contains the regulator e, as in cr = (1 — Co- e)/2, and 
that the limit Co- — > after integration is non-singular, since the exponents 
a and 7 in Eq. ( PP] ) are never half-integers. 

At two loops, the situation is similar, as may be seen from the expressions in 
Appendix 0. If Ci is the number of covariant denominators containing q, but 
not k, and C2 is the number of covariant denominators containing k, but not 
q, and if, furthermore, hi and 62 are the powers of go and ko in the numerator, 
respectively, and is defined as in Eq. (^S]), then a non- integer dimension of 
the energy components go, ko is mandatory in the following cases (i G {1, 2}): 

Ci — bi/2 + is < for ci and C2 7^ , 

Ci — bi/2 < for = 1 and Ci or C2 = . (21) 

As an example, consider the case is = 0, Ai = l,pi = l,Aj = for z > 1, 
Pj = for j > 1, and all ai,Pu = 1 in Eqs. (^2D - (^). (More general 
propagator powers do not change the arguments.) In that case, we have 

M4=(^l' ; (DetA^4)"^ = (xiX2)-^. 

Analogous to the one-loop case, additional powers of go or /cq in the numerator 
lead to more negative powers of xi or X2, as can be seen from Eqs. (0). For 
example, two powers of go or ko in the numerator lead to terms such as 

{DetM^y-^JillPli'^] ^ {DetM^y^M^^y^'^ Ji[l] 

^ 6,, f'dx.x-'-" , z,jG {1,2}. (22) 
Jo 

^By complex we mean that a has to be understood as a parameter in the complex plane 

in order to define the analytic continuation of the corresponding F-function. 
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Although higher powers of go or in the numerator will lead to even more 
negative powers of Xj, such as the parameter 7 can never be half- 

integer. As a result, the parameter integral is well defined for complex values 
of (T, and the limit a — > 1/2 after integration exists due to the analytic 
properties of the F-function. 

These arguments and the conclusion that 7 will never be half-integer, can 
be generalized to L loops by using the same formalism as in Appendix 0, 
with the objects UyAf^, . . . and l,Af, . . . being defined in L ■ 2a and L ■ 2u 
dimensional space, respectively. Hence, it follows that a troublesome F- 
function of the form F(±fc (1/2 — 0")) never arises. 

The arbitrariness in the choice of is not problematic, since it is just the 
noncovariant version of the arbitrariness in the definition of the subtraction 
scheme in the context of renormalization. To illustrate this point, consider a 
simple one- loop integral T{p) in split dimensional regularization: 



T{p) 



{2tiY {p-qfq^ 



(p^) 2-^F(2 - -) Beta(- - 1, - - 1 - a) . (23) 



(271)"'"^ ' ' 2' '2 '2 



Using n = 4 — 2cne, together with Eq. (|T3), and expanding in e, we arrive 

at 

T(p) = -2-^(p2)--|— + 4-21n(2)-^ + ln(47r)-74 • (24) 

The factor of l/c„ in the pole part is irrelevant since it will be removed 
by renormalization. As the terms ^(4-7?) and 7^;, in the subleading part of 
the e-expansion, are always present as an artifact of standard dimensional 
regularization, the additional term Co-/c„ is just another aspect of the arbi- 
trariness inherent in renormalization, caused by the breaking of covariance. 
The freedom in the definition of the finite part merely reflects the fact that 
renormalization schemes need to be specified. The procedure of dealing with 
the parameter during renormalization (i.e., whether or not it is subtracted 
together with the terms ln(47r) — 7^;) defines a specific "noncovariant renor- 
malization scheme" in the context of split dimensional regularization. 
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Of course, the easiest option would be to set to zero after having evalu- 
ated the integral, since, as explained above, the operation of setting Co- = 
ajter integration is well defined, and is the natural choice as it corresponds 
to (i(go) = 1 and di^q) = n — 1 in standard dimensional regularization. Nev- 
ertheless, it may be advantageous to keep the parameter useful check 
when computing gauge invariant quantities, since c^, obviously has to cancel 
out in any physical quantity. 

3 Leading divergence of the two-loop quark 
self-energy 

There are two methods of determining the overall divergence of the two-loop 
quark self-energy, namely by 

1. direct calculation of the two-loop integrals corresponding to the graphs 
shown in Fig. |l|; 

2. computing the counterterm diagrams shown in Fig. ^ by first extracting 
the UV-divergent parts of the corresponding one-loop subgraphs, and 
then calculating the divergence of the overall diagram containing these 
one-loop insertions. 

We have used both of these methods in order to check the consistency of our 
results. Note that method (1) requires us to keep the full e-dependence from 
the first integration if the momentum integrations are performed sequentially 
(this procedure is sometimes called the "nested method" ||19|)- Since it may 
not always be possible to maintain, without expansion, the full e-dependence, 
the majority of our two-loop integrals has been computed by integrating 
over both loop momenta simultaneously (also referred to as the "matrix 
method" [|r^), as outlined in Appendix]^. 

It is worth emphasizing that the (spurious) infrared divergences are more 
severe in the Coulomb gauge than they are, for example, in the Feynman 
gauge, the severity being caused by the noncovariant factor 1/g^ in the gluon 
propagator. A partial list of the noncovariant integrals used in our calculation 
is given in the Tables in Appendix ^ and |B|. The number of integrals per 
graph obtained by method (1) is typically of the order of a few hundred. The 
algebra has been performed by using the symbolic manipulation program 
FORM [13 . 
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Application of method (2) requires insertion of the following one-loop coun- 
terterms ^ ^ (cf. Fig. |): 



(25) 





= ^Cf-- - p , 
47r e 


-pa 


= nrw + rf)), 


rib) 


= i[Cf- —)—-[ 
2 Ait £ 


r(c) 






47r e L 



4 gn N 8 

3 ra^ 3 ra"' 



n 



F,ab 



47r e 



4 



{Nl - l)/(2iV, 



Tf = TR- Hf 



Note that S(p) and satisfy the Ward identity 

(p'_p).(r(f) + rW) = s(p)-s(p') 



(26) 



and that the divergent part of the ghost contributions, discussed in ref . ||T6[ , 
is zero. 
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p ^ 

(a) 




(be) 



(a) 



-4- 




+ 




(b) (c) 



(d) (e) 



Figure 2: Counterterms for the two-loop quark self-energy. The grey circles 
denote the UV divergent part of the corresponding one-loop insertions given 
below. 
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The two-loop corrections to the quark self-energy shown in Fig. |I| are given 
by the expressions in Eqs. (^9D to (|3TD , where we use the following short-hand 
notations: 

ci = + ie ; C3 = {p - ki)^ + ie , 
C2 = kl + ie ; C4 = {p - ^2)^ + ie , 

C5 = {k,~k2y + ie, (27) 

and 

Vr^^'^'igu q2, qs) = g^'^'ig^' - q^') + 9^'^'{q^' - q^ + 9^'''{q^' - qH , 

2 

d^^viq) = 9t,v + -^qi^qv - {qiJ^v + n^,qy) . (28) 
The "rainbow graph", diagram (a), is given by 

•c;(p-fci)C"(^2), (29) 

while the graphs with the Abelian and non-Abelian vertex subgraphs have 
the form 

(27r)^" 2 J C1C2C3C4C5 

•c;(p-fci)C"(^2), 



27r 2" 2 J 



(27r)2" 2 y C1C2C3C4C5 
■ C'^(fci) (^1 - C;(^2) V^37'(^i' -^2, A;2 - h) . 
The gluon self-energy insertion reads 



(27r)" 2 k^{k-qY 
■yiT^il^ ^ ~ ^) ^3'/''(-9> k,q-k) , 

such that 

sW(p) = ^T'^T'j d-k, ^^^^~2^^"^^' nr^(fco c;(fci) ■ 

(30) 
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The fermion loop insertion in diagram (e) is trivial (it does not contain a 
gluon propagator in the inner loop), and is given by 



n 



T\2 



ig' f rf-fc, n!:r (fci) C/(fci) • 



(27r)'^ 



C1C3 



(31) 



We would like to point out that S'^'^^ and E*^*^^ contain the double noncovariant 
denominator l/(/c^)^ stemming from the two terms of ^^"'^^(fci). In diagram 
(e), those double noncovariant denominators are reduced to single ones by 
corresponding terms in the numerator, whereas no such cancellation occurs 
in diagram (d). The integrals for graph (d) containing a factor l/{kf)'^ are 
those given in Table ^ for (3 = 2. 

After inserting the appropriate integrals and performing various crosschecks, 
as outlined above, we obtain the following results for the leading poles: 



Hp) 
Hp) 
Hp) 



2^ 



4 
27r 



e2 4 ' 



271 

27r 



e2 4 



41 
36 



9 



1 (I 



2txJ 



Note that the Ward identity ([26|) manifests itself in the equation 



(32) 



(33) 



We also recall that the one-loop gluon self-energy, n|^^"''(g), is not transverse 
in the Coulomb gauge, since it satisfies the Ward identity [|TB| 



q>'Il':f{q) + {q^g,, - q,q.)H-'^ = 



(34) 
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where H"''''^ is a non-vanishing ghost term. It should come as no surprise, 
therefore, that the final expression for T.^'^^ in Eqs. ( |32D still contains the 
vector n^. 

Finally, we should like to draw the reader's attention to the complete ab- 
sence of nonlocal factors in any of the leading-pole expressions, Eqs. (^2|). 
The absence of terms such as Po/p'^,p'^/p'^, etc., which implies the cancella- 
tion of all spurious infrared divergences, is certainly remarkable in view of 
the ubiquitous appearance of nonlocal terms at intermediate stages of the 
calculation. Of course, the trend was already set at the one-loop level, where 
the divergent parts of both the gluon self-energy and the quark self-energy 



were shown to be local [16, 17 



4 Conclusion 

In this paper we have tested the technique of split dimensional regularization 
by calculating the leading divergence of the two- loop quark self-energy 
in the Coulomb gauge. We find that the application of split dimensional 
regularization enables us to compute all two-loop integrals consistently, and 
that the final expressions for the leading poles of each graph are local, despite 
the presence of nonlocal terms at intermediate steps of the calculation. 
The leading pole parts of noncovariant integrals at one and two loops, used 
in the present calculation, can be found in Tables (Appendix ^, ^). We 
have also derived general expressions, in Feynman parameter space, for both 
one- and two-loop Coulomb-gauge integrals (Appendix |C|, 0). 
The latter expressions are useful in analyzing some general properties of the 
technique of split dimensional regularization. First of all, these expressions 
allow us to identify the class of integrals for which a complex regulator of 
the form cr = (1 — Co-e)/2 is mandatory in order that these integrals be at 
all well defined. Furthermore, it turned out that the leading pole of a tj^i- 
cal Coulomb-gauge integral, evaluated with split dimensional regularization, 
depends only on the sum of the regulating parameters, i.e. on a + uj = n/2, 
whereas the subleading pole generally contains the parameter c„. The lat- 
ter stems from the isolated appearance of the regulator a and is a direct 
consequence of breaking covariance. Finally, we have demonstrated that the 
limit Co- — i> after integration always leads to non-singular expressions, thus 
establishing the consistency of split dimensional regularization in general. 



19 



Acknowledgements 

We should like to thank Jimmy WiUiams for several discussions. This re- 
search was supported in parts by the Natural Sciences and Engineering Re- 
search Council of Canada under Grant No. A8063. 

A Divergent parts of one-loop integrals 

In the following tables we give the results in Minkowski space for the poles 
of the one-loop integrals entering in the calculation. All tables have to be 
read as follows: The result for the pole part of an integral 

ZTT^ J B 

is listed such that the denominator B is given in the first row and the cor- 
responding numerators A in the first column. All entries are implicitly mul- 
tiplied by r(e). Note that the poles which are of an infrared nature have a 
negative mass dimension (for example, a coefficient proportional to 
or where a > 0). The symbol "-" means that the corresponding 

integral docs not occm^ in our calculation and hence has not been compTited. 



B 

A 


q^{p — q)^q'^{p — qf 


q^{p — qY{k — qfq'^ 


{p — — qfq^ 


1 

Qo 


2po/{p^p^) 
2piPo/{p^P^) 

2prPj/{p'p'') 


Vip'k') 













-1 

1/3 5y 



Table 1: Pole terms of one- loop integrals in Minkowski space. 
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B 

A 


- q)^ 




a 


1 


1 + e 


2 + e 


1 


1 + e 


2 + e 


1 

Qo 

mo 

% 

1 

mj 


-2 
U 

/I 'r. /Q 



-2pV3 
-2pV3 
-16/15 + 
2/15 p2%. 


-1 
U 



-pV3 
-pV3 
-8/15piPj+ 
1/I5p2 5i,- 


-l/(2p2) 

U 

n 
U 



-1/2 
1/2 
1/6 


2/p2 




-1 
1 


l/p2 




-1/2 
1/2 
1/6 (5y 


[l + 4pV(3/)]// 
Po/p^ 
Pi/ (3p^) 


1/(2/) 
-1/(2/) 










q^'^ip — qY 


a 


1 


1 + e 


2 + e 


1 


1 + e 


2 + e 


1 

<lo 
Qi 

QiQo 

qI 

q 

QiQj 


2/p2 




-2 
-2 
-2/3 


l/p2 




-1 
-1 

-1/3 




1 

Po/2 
P^/2 
PiPo/3 
p2/3-pVl2 
p2/4 + pV3 


1/2 

Po/4 

Pi/4 

PiPol'o 
p2/6-pV24 

pV8 + pV6 
pV24 5ij+PiPj/6 


-l/(2p2) 




1/8 
-3/8 
-1/8% 



Table 2: Pole terms of one-loop integrals. All entries have to be multiplied by the overall factor r(e). 



B 

A 








a 
P 


1 

e 


1 

1 + e 


1 

2 


e 
2 


1 + e 
2 


e 
1 


1 

e 


1 

qo 
4 







-1 
-Po 

pV3-4p2/3 






-2 

-6po 
2p2 - I4p2 







[4 + 16pV(3p2)]// 
2po// 


1 






-1 


-2pi/3 
-2pVl5 




Table 3: Pole terms of one- loop integrals with non- integer denominator pow- 
ers. 
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B Leading poles of two-loop integrals 



The integrals listed are of the form 

where the denominator B is given in the first row and the corresponding 
numerators A in the first column. For the denominators, the short-hand 
notations defined in Eq. p7|), as well as fcs = p — fci, have been used. Note 
that fco = Po = P^ ■ AH entries are implicitly multiplied by 



B 

A 


C1C2 0304^2^/^3^ 


C2C'iC^kik^ 




CiC2C^k2k^ 


010204/^2^/^3^ 


1 

{pk2)kl 
{pk2)m 





1 








-4/p2 

-1 


-l/p' 



Pi 


-4/p2 


2p0 
2pl 




C'iC^k^k^ 




C\C^kc2 k^ 






1 

uo 

3 

Ph 


-2 

-2p° 





4 

4p0 

4/3(p2 + 2pl) 


2 



2p° 
-2/3 p 2 

2/3ip' + 2pl) 







Table 4: Leading pole terms of two-loop integrals with two noncovariant 
denominators. 
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B 

A 
/i 


CiC2C3cJ^2 


ClC2C3C5/C2^ 


CiC2C4C5^2^ 


C2C3CiC^k2 


1 

h.Oh.0 

{Vk2)kl 
{pks)kl 
{k2k3)kl 








-p72 





-p72 






1/4 

-p72 
p74 




1/4 

-p72 
p74 




C2C3C5k2 


CiC3C4/>;2 


^ ^ ^ ^2 
C3C4C5/t2 




1 

1.0 

"'2 

^3 
Pk2 

pks 


1 






pV2 


-2 
-2pO 

-2pV3-4p73 

_p2 


-1 

l/6p2 

-p73 - 2^2/3 
pV3 





Table 5: Leading pole terms of two-loop integrals with one noncovariant 
denominator. An overall factor 1/e^ is implicit. 

C General formula for one-loop integrals in 
the Coulomb gauge 

Here, we give the expression in Feynman parameter space for scalar one- 
loop integrals in the Coulomb gauge with an arbitrary number and arbitrary 
powers of propagators. We also outline how to obtain expressions in Feynman 
parameter space for the tensor integrals. 

Consider a scalar integral with r covariant and m noncovariant denominators 
in n dimensions: 

f (Pq 1 
Ir+m{aj,(3i) = y-^nn^^oT^' (35) 



24 



Going to Euclidean space and introducing Feynman parameters in the usual 
way, we find that the integral Ir+mi<^j, A) is given by 

'n „ r m 



I- M„ r m -1 



1 



a, — 1 



r=i j^i f^i J 

where 

r m 

r r r m 

V = g 2 ^ ^ Xj- g| - 2^4 ^ Xj P4J - 2q Xj pj + ^yiSi) 

3=1 j=l j=l 1=1 

r m 

j=i 1=1 

In order to obtain a formula for general tensor integrals, we first define 



-f n n ' (36) 



where is an arbitrary Lorcntz vector in Euclidean space, and then derive 
the tensor integrals by differentiation with respect to a: 



J ^7^2 ts, , , 



d d 



{-ir^---^Ir+m{a,a,,Pi) . (37) 



Carrying out the momentum integration in 7r+m(Q, ctj, Pi), we obtain 

= nu rfe) nE: r(A) H ft 
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Xj 



m r m „qq 

n dm vt" 5(1 - E - E m) \ dz ^^('^./^)— 1 
i=\ 1=1 i=\ •'^ 



\2 r 



exp I - 2;[ 2^ P4^^. — — + XjPj 

j=l ^J=l -^J j=l 



m 

\2 



+ E ^"'^ - (E ^jPj + 12yi ^if 
1=1 j=i 1=1 

■ exp{-/(a4)} • exp{-g{a)} , (38) 

/M = I Q4E a;iP4,j - I ' 

Z^j=i y j=i z ^ y 

/ r m \ la^ 

g{a) = a E ^i^i + E ^0 ~ ~T ■ ^^^^ 

Note that 7(04) and 51(a) depend on z, leading to more severe UV divergences 
for higher rank tensor integrals. 

For a scalar integral, we can immediately set a = 0, carry out the integration 
and go back to Minkowski space to get 

7. (a d,) = (-i)E,". r(A(a,/j)-a-a;) f fx . ^-.-^ ( y ^ . 

m r m 

l[dyiy?'~'^i^-j:^j-j:yi) 

1=1 j=i 1=1 

[-l^Xj POJ + 'r ^ + E ^jPj 

j=l 2^3=1 -^J j=l 



m 



+E?/'«' -ii^^jPj + i^yi^i) ■ (41) 



;=i j=i 1=1 
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D General formula for two-loop Coulomb gauge 
integrals 

A general two-loop integral in the Coulomb gauge, with r + m + 1 covariant 
and a + b + 1 noncovariant denominators, is of the form 

r,m,a,b - J .^^^ _ ^^^^^ _ ^^2(^ _ ^.)2A, _ p^.)2p, ^^)2a,(^ _ ^J2^„ " 

(42) 

The general exponents Xi,pj,ai and Pu have been introduced in order to 
account for certain denominators such as {{q — p)"^) , occurring for example 
in the gluon self-energy correction. 

Going to Euclidean space and introducing Feynman parameters Xj for the 
covariant denominators and yi for the noncovariant denominators, where xq 
and I/O are associated with {q — k)"^ and (g — k)"^, respectively, we find that 



r r+m a a+b 

/ n^^^^"'"' n dx^x^^-'j{dyiyt^-' [J dyuy'^-^ 

■'^ i=0 j=r+\ l=Q u=a+l 

r+m a+h 

S{l-J2xj-Yl yi) / dz . Ji[l] , (43) 

3=0 1=0 



where 



/n 1 
^ exp I - z[hMih + IMl- ^Nih - 2A^r+ 7^] | 

- ^^P^-^^> y'rf4.;^cxp{-^[/4W4-2A4/4]} 



2{a+Lo) 



j dHexp { - z[lMl- 2^1]} , (44) 



y —Xq Xo + Ej=H-l 
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M 



Ej=o + Et=o yi -xo - yo 

-xo - yo xo + yo + Y.'jtT+i Xj + Y.t=a+i yi 



'}2j=r+l Xj Pij 



= I t^i 



E"j=iXjPj + Et=iyisi 

Si=r+1 Xj Pj + Y.l=a+1 yi '''I 



■'3 

r+m a+b 

n = Y,Xjp] + Y.yi^i 

j=i 1=1 



2 



r r+m a a+b 

^{r,m,a,b) = XI -^^ + H Pj + H + H ' Aq = ao = 1 • 

i=l j=r+l 1=1 u=a+l 

Observe that for the generic case Aj = pj = ai = Pu = ^, (,{r,m,a,b) 
r + m + a + b. Integration over I in (Hf) leads to 



Ji[l] = (DetA^4)"MDetA<)"2 

■exp{-z[n-Af4M^^Af4-AfM-^Ar\}, (45) 



where 



r+m r r+m „ 

"1 Z(7 



2lo 



BetMi = [xo- J2 + J2 ' (46) 

j=l j=l j=r+l 

r+m a+b r a r+m a+b 

DetM = [ixo + yo)iY.^j + T.yi) + iT.^j + T.yi)iT. E yd 

j=l 1=1 j=l 1=1 j=r+l l=a+l 

By differentiating the integral ( PD repeatedly with respect to A/4 and/or Af, 
we may easily derive the appropriate expression for any tensor integral. 
Denoting the integrals with non-trivial numerators by Ji [l^^] , Ji [Z^*-'] , Ji [14^^^], 
etc., where l^^ = q^jf^ = k^,!^^^ = q,l^'^'' = k, and applying result P5|), we 
get: 

-^Ji[l] = 2zl?Ml] = 2z Jj/f] = 2z (Af.M^'Y^ Ml] . 
The procedure for other integrals is similar, leading to expressions such as 



Ml^^] = (AfM-'fMi]., 
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M^P^] = { {AfM-f {AfM-f + 1 {M-f ' }M1] , 

+l{M,f\^M-f'}Ml], (47) 



To cope with integrals in which either the factor {q — fc)^, or (g — k^, is 
absent in the basic integral (^), we proceed by first introducing the indicator 
functions 

r 1 if (g — fc)^ is present . r 1 if (g — fc)^ is present 

I if (g - A;)^ is not present ' ^"^ I if (g - A;)^ is not present 

(48) 



and 



6{xo) for = . . _ f 6{yo) for zc = 



md(^^) = |^ forzH = l ' ^^d(zc) = |^ 



for 2^ = 1 ' I 1 for ic = 1 ' 

such that xo is set to zero if {q — is absent, and yo is set to zero if (g — k^ 
is absent. After integration over z, the most general formula for the scalar 
(basic) integral (|i3|) is then given by 



7 = ( i)y.^^+y. p^+^B Tj^jr, m, a, b)+iB + ic- 2a -2uj) 

^ ' ^ ^ nLin-t";inr=ire±Uir(A.)r(p,)r(«or(/3.) 

„l r r+m a a+b 

/ Wdx^x":^-^ n dx,x'/-'l[dyiyr'-' J] ^?/«?/^"' 

i=0 j=r+l 1=0 u=a+l 

r+m a+6 

5(1 - -T.yi)- ind(^ij) ■ ind(zc) ■ (DetA^4)'^ (Det^)"^ 

j=0 1=0 

n - UiMl^U^ - UM-^U] . (49) 

The result in Minkowski space is obtained by making the replacements ip4, 
Po and pI; —p^. 
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Although the expressions given above look somewhat complicated, since they 
are intended to represent all types of integrals that may occur in a two-loop 
calculation, they can be easily implemented into an algebraic manipulation 
program in order to yield the parameter representations of the various inte- 
grals. 
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